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INTRODUCTION 


In the design of spacecraft structures, fine tuning the 
structure to achieve minimum weight with natural frequency 
constraints is a time consuming process. In this paper, a novel 
implementation of the method of optimality criterion (OC) is 
developed. In this new implementation of OC, the free vibration 
analysis results are used to compute the eigenvalue sensitivity 
data required for the formulation. Specifically, the modal 
elemental strain and kinetic energies are used. Additionally, 
normalized design parameters are introduced as a second level 
linking that allows design variables of different values to be 
linked together. With the use of this novel formulation, synthesis 
of structures with natural frequency constraint can be carried 
out manually using modal analysis results. Design examples are 
presented to illustrate this novel implementation of the optimality 
criterion method. 
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PROBLEM STATEMENT 


The optimal design problem to be solved is determination 
of the values of design variables such that the structure weight 
is minimized while maintaining a specified fundamental natural 
frequency of the system. The design variables are sizing of 
structural members, e.g. cross-sectional areas of truss elements; 
area moment of inertia of beam elements; and thickness of plate 
elements. Bounds on design variables are also considered. 
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OPTIMALITY CRITERION 


The optimal design problem defined in the previous section 
can be solved by mathematical programming techniques. To derive 
a simpler approach, we will treat the frequency constraint 
defined in Eq. (2) as equality constraint. Additionally, the 
side constraints will be ignored for the time being. With these 
simplifications, the set of optimality criteria can be derived 
among the design variables by the Lagrange multiplier technique. 
The optimality criterion can be interpreted as: 

At the optimal design, the ratio of the eigenvalue 
sensitivity to weight sensitivity is a constant for 
all design variables. 


FROM LAGRANGIAN : 


L = W - JX ( A, - A , J ) 

THE OPTIMALITY CRITERION: 
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EQUATION (5) IS THE OPTIMALITY CRITERION 
AN OPTIMAL DESIGN MUST SATISFY. 
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BASIC REDESIGN ALGORITHM 


Following the optimality criterion method suggested by 
Khot 11], linear recurrence relations can be developed based 
on Eq. (5). The Lagrange multiplier is computed by requiring 
that the updated design satisfy the frequency constraint. The 
basic redesign algorithm is summarized in Eqs. (6) and (7). 
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APPROXIMATE EIGENVALUE SENTIVITY ANALYSIS 


In the redesign algorithm, we need to know the derivatives 
of weight and eigenvalue with respect to the design variables. 
While the weight sensitivity is simple to calculate, the compu- 
tation of eigenvalue sensitivity could be quite involved because 
of the need to know derivatives of element stiffness and mass 
matrices with respect to design variables. In this paper, we 
adopt an approximate approach for computing eigenvalue sensi- 
tivity which use elemental strain and kinetic energy in the 
vibration mode [2]. 


EIGENVALUE SENSITIVITY 


GENERAL EQUATION: 
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M| W ^ ax. 

SIMPLIFIED EQUATION [2]: 
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ASSUMPTIONS: 
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SPECIAL CASE FOR: 

1. TRUSS ELEMENTS 

2. SYSTEM MASS MATRIX DOMINATED BY 
NON-STRUCTURE MASS 
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A NOVEL IMPLEMENTATION OF THE BASIC REDESIGN ALGORITHM 


The redesign algorithm given by Eqs. (6) and (7) can be 
implemented easily for truss structure. For example, if one 
uses MS C/NAS TRAN [3] for structural analysis, the strain energy 
and strain energy density can be obtained together with the 
modal analysis results. Using these data and assuming that a 
majority of the system weight is contributed from non-structural 
mass, the redesign algorithm can be implemented using the 
following procedures. 


(1) PERFORM MODAL ANALYSIS WITH STRAIN ENERGY AND 
STRAIN ENERGY DENSITY CALCULATIONS. 

(2) COMPUTE C : = USING 

** 3 A A, 

C x = I (3 e (ESE) e /CESED s .X e ) 


(3) COMPUTE G^ = USING 

<7 A A. 

Qa = 2 I (ESEe/XeM.) 


NOTE THAT THE SUBSCRIPT e REFERS TO ELEMENT NUMBER 
AND THE SUMMATIONS IN EQS. (13) AND (14) ARE OVER 
ALL THE ELEMENTS THAT ARE ASSIGNED AS DESIGN VARIABLE 
Xi. 

(4) COMPUTE LAGRANGE MULTIPLIER USING EQ . (7). 

(5) UPDATE THE DESIGN VARIABLES USING EQ . (6). 
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DESIGN EXAMPLE 


The design procedure described in this paper has been applied 
to a truss structure shown below. The objective is to find the 
minimum weight design of the truss structure to maintain a speci- 
fied fundamental natural frequency of the system. Starting from a 
uniform truss structure that satisfied the 5.7 Hz constraint on 
the fundamental natural frequency, the design is manually opti- 
mized by typical trade studies. The optimality criterion algo- 
rithm is then applied to this manually optimized structure to 
obtain an additional 25 pound saving in the structural weight. 

The comparison of the truss structures weights is shown in the 
table below. 


DESIGN CASE 

TRUSS 

STRUCTURE 

WEIGHT 

f 

UNIFORM 
TRUSS SIZE 

2102.0 

lbs 

5.7 

Hz 

MANUALLY 
OPTIMIZED 
TRUSS SIZE 

1014.0 

lbs 

5.7 

Hz 

OPTIMIZED 
TRUSS SIZE 
USING METHOD 
OF OC 

989.0 

lbs 

5.7 

Hz 
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CONCLUDING REMARKS 


A method of optimality criterion was shown to be a powerful 
tool for minimum weight design of structures with constraint on 
fundamental natural frequency. With the new method of implementation 
presented in this paper, the design procedure can be carried out 
by simple calculations. The effectiveness of this approach has 
been demonstrated by a truss structure. This method can be extended 
to other types of structure elements using eigenvalue sensitivity 
formulation in Ref. [2J. 
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SYMBOLS AND ABBREVIATIONS 
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aw/jy; FOR current design 

STRAIN ENERGY 
STRAIN ENERGY DENSITY 
FIRST NATURAL FREQUENCY 
DESIRED FIRST NATURAL FREQUENCY 

aXi/ax; for current design 

GLOBAL STIFFNESS MATRIX 
GLOBAL MASS MATRIX 

GENERALIZED MASS OF THE FIRST MODE 
SPACE OF DESIGN VARIABLES 

TOTAL KINETIC ENERGY OF ELEMENTS ASSOCIATED WITH 
DESIGN VARIABLE i FOR MODE 1 

TOTAL STRAIN ENERGY OF ELEMENTS ASSOCIATED WITH 

DESIGN VARIABLE i FOR MODE 1 

STRUCTURE WEIGHT 

DESIGN VARIABLE VECTOR 

LOWER BOUNDS OF X 

UPPER BOUNDS OF X 

CURRENT DESIGN VARIABLE 

CURRENT DESIGN 

UPDATED DESIGN 

EIGENVALUE 

RELAXATION FACTOR 

WEIGHT DENSITY 
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